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ABSTRACT 


A  computer  program  has  been  developed  to  compute 
initial  satellite  orbits  from  two  or  more  independent 
observations  of  satellite  position.  The  technique  combines 
other  observations  with  the  initial  elements  by  means  of  a 
differential  correction  scheme  to  obtain  a  corrected 
initial  orbit.  Extensive  experimentation  indicates  that 
appropriate  selection  of  observations  to  be  used  in  the 
initial  computations  will  yield  accurate  initial  orbits. 
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SECTION  1 


INTRODUCTION 


The  Initial  Orbit  Radar  Fix  Pro 
( I0RF1  is  coded  for  the  Phiico  2000  compu 
of  the  Semi-Automatic  Programming  System 
Springs,  Colorado.  It  functions  in  conju 
and  receives  its  input  from  ihe  Schedule 


gram  with  Differential  Correction 
ter  and  is  designed  to  be  a  part 
at  the  SPADATS  CenLer  in  Colorado 
nction  with  the  Executive  Program 
and  SEAI  Tapes. 


The  purpose  ut  the  program  is  to  compute  a  set  of  initial  elements, 
for  a  satellite  orbit,  from  two  or  three  position  fixes.  These  position  fixes 
may  be  isolated  radar  strikes  which  come  from  different  stations  or  occur 
during  different  revolutions.  The  method  used,  is  to  make  a  first  approxi¬ 
mation  of  the  semi-latus  rectum,  p,  and  iterate  until  the  computed  time 
between  the  two  observations  is  within  some  specified  tolerance  of  the 
observed  time  between  the  two  observations  (p- i tera t ion) . 


The  option  to  differentially  correct  the  initial  orbit  is  avail¬ 
able.  The  differential  correction  technique  used  is  the  same  as  that  used 
in  the  SGPDC  Program.  The  inclination  may  be  corrected,  but  the  drag  co¬ 
efficient  is  never  corrected.  Regardless  or  whether  or  not  the  option  Lo 
correct  the  orbit  is  chosen,  the  seven  standard  element  cards  are  always 
punched . 


SECTION  2 

FORMULATION 


2 . I  COMPUTATION  OF  r  FROM  A  POSITION  FIX 


Given  the  quantities: 

Observation  day  number.  (in  days  since  1950.0); 

Minutes  from  Greenwich  midnight  of  observation  day  to 
observation  time; 

Observed  azimuth,  A,  and  elevation,  h;  or  topocentric 
right  ascension.  Of  ,  and  declination,  6  ; 

Observed  slant  range,  P; 

Latitude,  0  ,  and  east  longitude,  X  E,  of  the  observing 
station; 

-(C+H)  cos  <£»  and  Z  =  -(S+H)  sin  v  for  the  observing 
station; 

The  Greenwich  sidereal  time  at  the  start  of  the  obser¬ 
vation  year,  >;grQ; 


obtained 
by  0BSGET 
routine 


obtained  by 
TLC  routine 
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The  following  procedure  is  used  to  compute  the  geocentric  position 

vector,  jr: 


a  Compute  the  geocentric  station  vector,  R. 

0  gr  (deg)  =  ®8r0  +  (  ® -  360°)  D  +  0F 

Where  0gr  is  Greenwich  sidereal  time  at  observation  time,  9  is  the 
rotation  rate  of  the  earth  in  deg/solar  day,  and  D  and  F  are  the  observation 
time  in  days  and  fractions  of  a  day,  respectively,  since  the  start  of  the 
observation  year. 

9  = 6  +  X_  where  9  is  the  observation  sidereal  time, 

gr  E 

X  =  [  -  (C+H)  cos  0  ]  cos  9 

R  Y  =  [  -  (C+H)  cos  0  ]  sin  9 

Z  =  -  (S+H)  sin  0 


b.  Compute  the  topocentric  position  vector,  p 

L  =  cos  9  (cos  0  sin  h  -  cos  A  cos  h  sin0)  -  sin  Scos  h  sin  A  ^ 


=  sin  9  (cos  0  sin  h  -  cos  A  cos  h  sin  0)  +  cos  9cos  h  sin  A 

L  =  sin  0  sin  h  +  cos  A  cos  h  cos  0 
z 


If  A,  h 
f’  are 
observed 


L  =  cos  6  cos  a 
x  / 

L  =  cos  6  sin  a 

y 


L  =  sin  6 
z 


J 


if  a  ,  6 
are  observed 


£=  PL 

3)  Compute  _r  and  r 
£  =  £-  R 


=  sT 


v/  2  ,  2  2 
=  V  x  +  y  +  z 


2.2  PRELIMINARY  COMPUTATIONS,  INDEPENDENT  OF  p 

Given  two  position  fixes,  r^,  and  £2,  the  following  computations  are 
made.  These  calculations  are  independent  of  the  selection  of  the  semi-latus 
rectum,  p. 
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Compute  the  Orientation  Vectors  and  Angles. 


U.  -=i  1-1.2 

—  l  r . 

l 


c  = 

cos  (v^  - 

V 

=  ”i  ' 

—2 

s  = 

sin  (v2  - 

V 

-  ±  V 

1  -  c2 

w  = 

Si  X  —2 

(Wz 

=  cos 

i) 

S 

•H 

>1 

=  W  x  U. 

—  —l 

i  = 

1,  2 

(both  signs  are  tried) 


sin 

u/T 

-  w2 

Z 

i  = 

tan  ^ 

sin  i 
cos  i 

0  <  i  <  7T 

a  = 

tan  ^ 

W 

X 

o  <  n  <  2tt 

l)  Calculate  Range  of  p's  and  Maximum  Number  of  Revolutions  Possible 

An  expression,  which  describes  the  curve  on  which  the  solution  for 
the  semi-latus  rectum  must  fall,  can  now  be  derived  using  the  following 

expressions ; 


e  cos  v. 

l 

e  sin  v^ 


e  sin  V£ 


p/ri  -  1  i  =  1 ,  2 

C  e  cos  v  -  e  cos  v 2 
S 

e  cos  v^  -  C  e  cos  v 2 
S 
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These  equations  may  be  substituted  into  the  expression 


to  yield 


2  2  2 
e  =  (e  sin  v. )  +  (e  cos  v  ) 


2  2 

e  =  p  +  a2  P  +  a3 


where 


a  =  _i—  +  — 

1  rx2  S2 


a2  =  2 


(*■*)’ 


-  7"  +  \  (c-i) 
ri  s2 


and 


a  =  1  +  ^  (C- 1 ) 2 

S 

2  2 

The  equation  for  e  describes  a  parabola  on  an  e  ,  p,  plot.  To 
determine  the  limits  of  p,  the  quadratic  may  be  solved  for  the  arbitrary 
limiting  values  of  e2  =  0.95,  (e  =  1  gives  singularity  points.) 


P  = 


-  -  4a1  (a3  -  0.95) 


2a 


1 


Therefore,  for  any  e  <  0,95  the  roots  of  the  above  equations  must  bound  the 

actual  value  of  p,  i  e,. 


Plmin  -  P  —  Plmax 


where 


D 

lmax 


a9  +y  a  2~  4a  x  (a3  -  0.95) 


2a, 


-a2  -  \[a2*  -  4ai(a3  -  0.95) 

2T  ‘ 


P  = 

lmin 
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The  subscript  1  is  used  to  differentiate  these  p's  from  the  ones  which  follow. 
Similarly : 

P2min  —  P  —  P2max 

where  p  .  is  derived  from  Kepler's  second  law 
'zmin 

P  ~  (ke(t2"tl))  =  P2min 

and 

p  =  r(l  +  e  cos  v)  <  2(the  smaller  of  ^  or  r 2)  =  P2max 

There  is  also  the  additional  requirement  that  p  be  always  greater  than  unity. 
Thus  the  maximum  allowable  p,  max  p,  is  taken  as  the  smallest  value  of  P^max 

and  p  Similarly  the  minimum  allowable  p,  min  p,  is  taken  as  the 

2max 

largest  value  of  1,  Plmin,  and  p2mln 

The  maximum  number  of  revolutions  possible  during  the  time  interval 
(t2-t1)  can  be  calculated  by  the  following  equation: 

C2  '  C1 
max  N  <  — - - 

co 

where  P  is  the  period  of  a  surface  circular  satellite, 
co 


2  3  p- ITERATION 

The  interval  for  p  (min  p  <  p  <  max  p)  is  divided  into  m  parts, 
where  m  is  a  constant  and  is  part  of  the  program  input.  These  divisions 
may  be  called  A  p,  i  e. 

^  _  max  p  -  min  p 

P  m 
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There  is  now  a  family  of  p's;  min  p,  min  p  +  ^P.  tnin  p  +  2  Ap  .  .  .  , 

max  p  The  p- iteration  is  now  commenced  by  calculating  a  functio  , 
f  (min  p)  for  all  the  integer  values  of  N  less  than  max  N.  (T 
derivation  of  fM(p)  will  be  shown  below)  This  same  procedure  is  then 
plrformed  for  pN  min  +  A  p,  etc.  .  g]L.  If  there  is  a  change  in  the  sign 
of  f^p)  for  successive  values  of  p  at  some  N,  say  between  pL  and  p2 
for  N*,  a  root  is  indicated. 

f  /Pj  *  P_2  |  ig  computea  and  its  sign  examined  to  determine  on  which 


side  of 


Pi  -  P2 


the  root  is  If  the  root  is  less  than 


Pi  +  P2 


Pi  +  P2 


,  then 


£N*(3Pl  4  ~  ^  iS  computed;  if  the  r00t  iS  greater  than  2  £  '  then 
f  ^  ^Pl  +  3P2)  is  computed.  This  halving  process  is  continued  until  a 
p  is  found  for  which  \f^p)  I  <  €3,  (  £3  is  an  inPut  t0  the  Pr°gram) ' 

Given  p  and  N,  fN(p)  is  computed  as  follows: 
a  Compute  the  Shape  and  Size  of  the  Orbit  and  the  Period: 


e  cos 

V  =  -2-  -1  i  -  1,2 

1  r . 

1 

e  sin 

Cecosv^-o 

cos 

V1  "  s 

e  sin 

e  cos  v ^  -  C  e 

cos 

V2  S 

where 

C  =  cos  (.y 2  ~  vp  1 

•h  . 

2 

e  = 

2 

(e  cos  v^)  4-  (e  sm 

vL)2  = 

a  = 

p/(l  -  e2) 

n  = 

,  /  3/2 
k  /a 

e 

P  = 

27T 

n 

.2  .  ..  \ 2 
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b ,  Compute  the  Mean  Longitude,  Li,  i  -  1,  2: 


( 1  +  cos  i )  cos  SL  ^ 
(1  +  cos  i)  sin  i ^ 
( 1  +  cos  i )  cos  i ^ 


=  u  .  - 

V  .  \ 

yi 

XI 

=  u  .  + 

V  . 

XI 

yi 

=  -  U  . 

-  v  . 

yi 

XI 

=  u  .  - 

V  . 

XI 

yi 

if  W  >0 
z 


if  W  <0 
z 


SL.  =  tan 

l 


-1  (  Sin  jl  \f  1  +  cos  1  \ 

I  cos  FT  /(  1  +  cos  i  J 


JL 


s  in  ( E  -  v )  . 
r .  i 

l 


(e  sin  v^)  (e  cos  v^) 
1  +  sTl  - 

2 


+  e 


(  e  sin  vi> 

-2-  cos  (E  -  v)  .  =  1  ~ 
rt  1  1  +  \I  1  -  e 


- .  ^  +  e  cos  vt 

+  V 


(E  -  v)i  =  tan 


-  1 


(p/r^  sin  (E  - 
(p/,ri)COS  ’  V’^ 


J 


e  sin  Ei  =  y  ^ 1  '  (e  sin  v  i\ 

L  =  I.  +  (E  -  v),  -  e  sin  E 
li  i  1 


c.  Compute  i  (p) : 

- - —  L2  -  Lx 


A  t 


com 


0  <  L2  -  Lx  <  2  7T 


f  (p)  =  At 
NVF'  com 


At  ,  +  P  •  N 

obs 


sin  v. 

l 
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2.4  COMPUTATION  OF  OUTPUT  QUANTITIES 

The  following  quantities  are  computed,  in  addition  to  quantities 
already  computed,  for  inclusion  in  the  output  once  a  solution  is  found. 

a.  Compute  a^,  ayN>  and  0)  : 


V  ,  e  cos  v,  +  U  ,  e  sin  v 
z  1 _ 1 _ zl _ 

axN  ~  sin  i 


U  ,  e  cos  v,  -  V  ,  e  sin  v 
zl _ 1 _ zl _ 

ayN  ~  sin  i 


(JO 


tan 


-1 


a 


xN 


b.  Compute  the  time  of  nodal  crossing,  TM,  closest  to,  but  after  t1 : 


(!  +  axN)  e  cos  En  «  +  e 

(!  +  axN)  e  sin  E^  =  -  ^ 1  -  e2  ayN 


En  =  tan 


-  1 


( 1  +  a  )  e  s in  E 
_ _ xN _ N 

(1  +  a  )  e  cos  E 
v  xN  N 


e  sin  E 


4 


~  ? 

1  -  e  a 


N 


1  +  a 


xN 


£N 


”n  en 


E„  -  e  sin  E 


N 
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I 


M,  =  L.  -  03  -  n  if  w  >0 

11  z 

M,  =  L,  -  CO  +  a  if  W  <0 
11  z 

“n  '  Mi 

tn  "  ci  +  - - - .  where  0  <  -  Mj^  <  2  n 

c.  Compute,  .q_,  h,  and  c": 


q  =  a  (1-e) 


h  =  ^^p1  w 
c"  (1/min) 


-  360  n2  c 
_ o 

2 


where  c  is  in  days 
o  J 


2.5  THIRD  VECTOR  CHECK 

If  a  third  position  fix  is  available,  _Tj ,  it  is  used  to  check  each 
possible  solution.  The  position  on  the  possible  solution  orbit,  _r^ ,  at  the 
time  of  the  3rd  position  fix  is  computed  using  the  Simplified  General  Pertur¬ 
bations  technique*,  with  t^  as  epoch.  The  difference  between  the  computed 

and  the  observed  position  is  determined,!  Ar  |  =  * J  (rj  -  £3) -(£3  -  r_3) 

and  compared  with  e^,  (  is  an  input  to  the  program).  If  |  ArJ  <_  then 

the  orbit  is  assumed  to  be  a  solution  and  is  printed  with  all  other  solutions. 


*See  Section  4.49,  Aeronutronic  Publication  U-1691,  Semi-Automatic  Programming 
System,  Programming  Document. 


SECTION  3 


FLOW  DIAGRAMS 


On  the  following  pages  is  a  flow  chart  of  this  program. 

Diamonds  indicate  decisions.  Arrows  indicate  direction  of  logical 
flow.  If  a  solid  arrow  and  one  or  two  dotted  arrows  emerge  from  the  same 
box,  the  solid  arrow  indicates  the  normal  flow  of  logic  while  the  dotted 
arrows  indicate  some  special  condition  which  is  described  next  to  the  arrow. 
Ovals  and  circles  indicate  entrances  and  exits.  The  oval  at  the  entrance 
to  a  subroutine  contains  within  it,  the  symbolic  name  of  the  subroutine. 
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Output  observation  number 
t,  £,  r  for  the  3  (or  2) 
observations.  Check 
m,  and  maximum  number  of 
solutions  for  0  and  set  in 
standard  values  for  any  that 
are  0 . 


Set  error  exits,  Set 
up  page  heading  for 
initial  orbit  output 


Enter  XSRCH  routine  3  times  to  input  3 
control  cards  from  C0NBUF  buffer  (sat¬ 
ellite  number,  element  set  number, 
satellite  name,  differential  correction 
flag,  N0>  c0,  maximum  number  of  solu¬ 
tions,  e^»  m>  Nf>  i  correction  flag) 


Got  • 

j  position 
S'vfixes  yet 


No  Station  i 

Enter  RBAR  routine  to  get  next  Coordinates  j 

observation  from  0BL0C  and  compute  or  p *  o 
_r  ,And  r  for  this  observatio'n . 


I  Observation 


''Are  there^ 
2  position 
fixes  / 


Write,  LESS  THAN  2 
P0SITI0N  FIXES,  THUS 
N0  0RBITS  CAN  BE 
DETERMINED 


Enter  PITER  routine  to 
compute  orbit(s)  that 
satisfy  the  2  or  3 
position  fixes 


Write:  N0  0RBITS 
F0UND  THAT  SAT¬ 
ISFY  THE  1ST  3 
(0R  2)  POSITION 
FIXES 


rap  up  output  > 
tape  and  return 
control  to  exec 
\utive  routine  ^ 
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SECTION  4 


OPERATING  INSTRUCTIONS  AND  INPUT-OUTPUT 


This  program  operates  under  the  control  of  the  Semi-Automatic 
Programming  System  (SPS)  Executive  Program  using  the  Schedule  Tape  and 
the  S  File  for  input.  A  "cold  start"  loads  the  Executive  Program  into 
memory  and  transfers  control  to  it.  The  Executive  Program  examines  the 
toggle  switches  and  determines  that  it  should  continue  operating  in  the 
Schedule  Tape  Mode.  The  data  cards  for  the  program  indicated  on  the 
SPSJ0B  card  are  processed,  placed  in  the  proper  buffers,  the  program  is 
loaded  into  memory,  and  control  is  transfered  to  the  program. 

For  the  details  of  SPS  operation  see  Aeronutronic  Publication 

U-1691. 

INPUT 


The  makeup  of  the  card  deck  that  is  put  on  the  schedule  tape 
is  as  follows:  (Below  is  an  example  containing  multiple  cases) 


1. 

J0B  card* 

2. 

Remarks  card* 

3. 

SPSJ0B  card 

4. 

Three  control  cards 

for 

case  1 

5. 

Any  number  of  observation  cards*  for 
(sensor  cards*  also  for  input  option 

case  1 
2) 

*See  U-1691  for  format. 
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6.  End  of  case  card* 

7.  Three  control  cards  for  case  2 

8.  Any  number  of  observation  cards* for  case  1 
(sensor  cards*  also  for  input  option  2) 

9.  End  of  case  card* 

10.  Three  control  cards  for  last  case 

11.  Any  number  of  observation  cards*  for  last  case 
(sensor  cards  also  for  input  option  2) 

12.  End  of  case  card* 

13.  End  of  job  card* 


Card  formats  are  as  follows: 


Card  Columns 

SPSJ0B  Card  1-8 

9-16 

17 


18 

19-72 

73-79 

80 


SPSJ0BAA  ■) 

I0RF/WA  f 

0  or  2  (input  option) 


means  blank 

0  means  sensors  are 
obtained  from  the  S 
file  while  2  means 
they  are  obtained 
from  sensor  cards 


0  (output  option) 

Any  identifying  comments  (alphanumeric) 
Reserved  for  expansion 
J  (card  type) 


*  Ibid. 
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Card 

Columns 

Contents 

Control  Card 

1 

C  1-3 

Satellite  number  (alphanumeric) 

4-6 

Element  set  number 

7-8 

Not  used 

Right 

A 

9-18 

Satellite  name  (alphanumeric) 

Adjusted  1 

1  19-22 

Not  used 

1  23-36 

Epoch  revolution  number 

37-78 

Not  used 

79 

DC  option  -  0  means:  do  not  differentially 

L 

correct,  4  0  means:  do  differentially 

correct . 

80 

P  (card  type) 

Control  Card 

2 

f  i*8 

9-22 

Not  used 

c  in  days  (if  not  known,  use  0) 

23-24 

M,  maximum  number  of  solutions  desired 

25-36 

£3  in  mins.,  the  magnitude  of  the  difference 
between  the  computed  and  the  observed  time 

Floating  ) 

between  the  two  positions,  must  be  less 

Point 

than  £3. 

37-48 

£4  in  kms.,  the  distance  of  the  3rd  position 
fix  from  the  computed  orbit  must  be  less 
than  £^. 

49-60 

m,  number  of  parts  to  divide  the  interval  between 
minimum  p  and  maximum  p  into. 

i  61-79 

Not  used 

^  80 

P  (card  type) 

Control  Card 

3 

1-22 

Not  used 

23-36 

Prediction  revolution  number  (see  DC  output) 

37-66 

Not  used 

67 

i  option  -  0  means  "correct  inclination," 

4  0  means  "don't  correct  inclination" 

68-79 

not  used 

80 

P  (card  type) 

All  columns  used  for  identification  (e.g.  Control  Card  1,  Columns  1-3)  may  be 
left  blank  if  the  information  is  unknown.  If  M,  £3,  £4  and/or  m  are  input  as  0 

or  all  blanks,  they  are  replaced  internally  by  1.0,  .001  mins.,  100  kms.,  and/or 
50.0  respectively.  The  decimal  point  may  be  punched  anywhere  within  the  floating 
point  field. 


c 
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Before  transferring  control  to  this  program,  the  Executive  Program 
does  the  following  operations: 

1.  The  3  control  cards  are  placed  in  C0NBUF  in  BCD  character 
format  (80  columns) 

2.  The  observation  cards  are  converted  and  placed  in  0BL0C. 

3.  The  sensor  cards  (input  option  2)  are  converted  and  placed 
in  SBL0C ,  or  (input  option  0)  the  sensors  associated  with 
the  observations  are  read  from  the  S  file  and  placed  in 
SBL0C 


OUTPUT 


At  the  top  of  each  page  the  comment  "INITIAL  0RBIT  DETERMINATI0N , 
P-ITERATI0N  METH0D",  the  date,  and  the  page  number  are  printed.  The  first 
comment  that  is  printed  after  the  page  heading  for  each  case  is  "THE 
0BS ERVATI0NS  T0  BE  USED  ARE  followed  by  a  listing  of  observation  numbers., 
time,  x,y,z,  and  r  for  the  first  three  (or  two)  observation  fixes.  These 
quantities  have  appropriate  headings  that  also  give  the  units  used. 


For  each  possible  solution,  i.e.,  each  pair  of  p  and  N  for  which 
I  fN(p)  I <£3 »  the  following  quantities  are  printed,  a,e,i,  0,  CD,  Tn,  Lt ,  axN, 
a  N ,  hx ,hv ,hz , tl ,N ,q ,p ,n,  and  c."  Following  this  is  a  line  containing  the 
results  of  the  third  vector  check.  If  |Ar  |<  £  this  line  reads  "MAGNITUDE 
0F  DELTA  R  BAR  =  , XXXXXXX  KMS ,  THEREF0RE  3RD  0BSERVATI0N  D0ES  FIT  THIS 

0RBIT"  If|Ar|>e^,  the  line  reads  "MAGNITUDE  0F  DELTA  R  BAR  =  .XXXXXXX 
+  y  KMS,  THEREF0RE  3RD  0BSERVATI0N  D0ES  N0T  FIT  THIS  0RBIT" .  See  Figure  1 
for  the  format  of  the  ou'tput. 


There  are  several  error  conditions  that  cause  error  comments  to  be 
printed.  They  are  described  below. 

1.  If,  during  the  input  of  the  three  control  cards  any  illegal 
characters  are  encountered,  "ERR0R  IN  INPUT  CARDS"  is  written. 

2.  If ,  at  any  time  during  the  program,  exponent  overflow  causes 
a  jump  to  location  0,  "EXP0NENT  0VERFL0W  FR0M  L0CATI0N  XXXXX" 
is  written. 

3.  If,  at  any  time  during  the  program  a  subroutine  error  causes 
a  jump  to  location  3,  "SUBR0UTINE  ERR0R  FR0M  L0CATI0N  XXXXX" 
is  written. 
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4.  If  there  are  less  than  two  position  fixes  available,  "LESS 
THAN  2  P0SITI0N  FIXES,  THUS  N0  0RBITS  CAN  BE  DETERMINED" 
is  written. 

5.  If  ti  =  t2,  "THE  TIMES  ARE  EQUAL  THEREF0RE  N0  0RBITS  CAN  BE 
DETERMINED  FROM  THE  1ST  TW0  P0SITI0N  FIXES"  followed  by 
error  note  11  is  written. 

6.  If  t2-tL  =  2  days,  "THE  TIME  BETWEEN  THE  TW0  P0SITI0N  FIXES 
IS  N0T  LESS  THAN  2  DAYS,  THUS  THIS  R0UTINE  WILL  ATTEMPT  N0 
0RBIT  DETERMINATIN'1  followed  by  error  note  11  is  written. 

7.  If  sin  (v2-v1)  =  0,  "THE  DIFFERENCE  BETWEEN  THE  TRUE  AN0MALIES 
IS  EQUAL  T0  0  0R  180  DEGREES  THEREF0RE  N0  0RBITS  CAN  BE  DETER¬ 
MINED  FR0M  THE  1ST  TW0  P0SITI0N  FIXES"  followed  by  error  note  11 
is  written. 

2 

8.  If  there  are  no  values  of  p  for  which  e  <  .95,  "ALL  VALUES 
OF  P  GIVE  E**2  >  .95  THE RE F0 RE  N0  0RBITS  CAN  BE  DETERMINED 
FR0M  THE  1ST  TW0  P0SITI0N  FIXES"  followed  by  error  note  11 
is  written. 

9.  If  max  p  < min  p,  "THE  MAXIMUM  ALL0WABLE  P  IS  STRICTLY  LESS 
THAN  THE  MINIMUM  ALL0WABLE  P  THEREF0RE  N0  0RBITS  CAN  BE 
DETERMINED  FROM  THE  1ST  TWO  POSITION  FIXES"  followed  by 
error  note  11  is  written. 

10.  If  in  the  subroutine  R00T,Ap  becomes  so  small  that  p  +  Ap 

=  p,  "R00T  INDICATED  BETWEEN  P  =  X.XXXXXXXX  AND  P  =  X.XXXXXXXX 
F0R  N  =  XX  BUT  C0ULDN"T  C0NVERGE.  LAST  P  TRIED  =  X.XXXXXXXX, 

F  SUB  N  OF  P  =  +  .XXXXXXXX  +  YY"  is  written. 

11.  If  at  the  end  of  the  initial  orbit  part  of  the  program  no 
solution  orbits  have  been  found,  "N0  0RBITS  F0UND  THAT  SATISFY 
THE  1ST  X  P0SITI0N  FIXES"  is  written. 

After  writing  any  of  error  notes  1  through  9  or  11  the  output  tape 
is  wrapped  up  and  control  is  returned  to  the  Executive  Program.  After  writing 
error  note  10  the  program  continues  normally  to  look  for  solutions. 
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If  differential  correction  is  not  desired,  the  program  punches 
the  7  standard  element  cards  at  the  time  of  nodal  crossing  closest  to, 
but  after  t^,  for  the  last  solution  orbit  found. 

If  differential  correction  is  desired  then  the  orbit  that  is 
corrected  is  the  very  last  solution  found.  The  inclination  may  be  corrected 
if  desired,  but  the  drag  coefficient,  c",  is  never  corrected.  During  the 
differential  correction  the  outputs  are:  all  the  residuals  in  the  observations, 
corrections  to  the  orbital  elements,  and  corrected  orbital  elements  for 
every  iteration.  The  7  standard  element  cards  are  punched  for  the  final 
elements,  both  at  the  time  of  nodal  crossing  closest  to  t^ ,  and  at  the  time 
of  nodal  crossing  into  the  prediction  revolution.  For  details  of  the 
above  output  method,  reference  should  be  made  to  the  description  of  the 
SGPDC  Program. 
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INITIAL  0RBIT  DETERMINATI0N.P-ITERATION  METH0D 


M 

M 


I 


DCJ 


I 


+  1+1  + I 


+  1+1  +  1 


+  1+1  + I 


+  1  +  1  +  1 


CO 

CO 

CO 


CO 

B 

£ 

£ 

M 

H 

M 


H 

CO 


SL 

Pm 


Pm 

S 

SI 


w 

Di 

D 

O 

M 

Pm 
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SECTION  5 


TEST  CASES 


Both  .cases  use  input  option  0. 

CASE  1  : 

INPUT: 

Control  Card  1 
Sat.  no.  =  163 
Element  set  no.  =  012 
Sat.  name  =61  SIGMA  1 
Epoch  rev.  no.  =  1288 
DC  option  =  1 
Control  Card  3 
Prediction  rev.  no.  =  1288 
i  option  =  0 


Control  Card  2 

c  =0.0  days 
o 

M  =  99 

e 2  ~  0-001  mins. 

£ .  =  100.0  kms . 

4 

m  =  2 
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Observation  Cards 


I 


Sta 

No. 

Time  of 

Obs. 

I  /I  Q  fr  i  /-\  y* 

A(deg)  or 
a  (hr , min, sec) 

Year 

Mo 

Day 

4 

Hr 

Min 

Sec 

niaefc, j  or 
6(deg) 

p(kms) 

0001 

1 

12 

2 

21 

27.0 

h  =  30.5 

A  =  10. 1 

4800.0 

0001 

1 

12 

4 

2 

11 

27.0 

h  =  1.7 

A  =  1.6 

7167.0 

0001 

1 

12 

4 

2 

51 

27.0 

6  =  40.8 

a  =  3him36?o 

6760.0 

0001 

1 

12 

4 

2 

41 

27.0 

h  =  38.0 

A  =  170.3 

4457.0 

0001 

1 

12 

4 

2 

46 

27.0 

h  =  20.2 

A  =  175.9 

5543.0 

OUTPUT : 

The  first  outputs  are  the  first  three  position  fixes 
Time 


YMMDDHH 

MMSS.SSS 

x(radii) 

y(radii) 

z(radii) 

r(radii) 

61120402 

2126.979 

0.403031251393 

0.440155155200 

1.40202623326 

1.52376123762 

61120402 

1126.981 

-0.0164250853704 

0.0267504054063 

1.52240997878 

1.52273356565 

61120402 

5126.993 

1.09301779186 

1.08379045734 

.0189203224144 

1.53936598217 

The  only  solution  printed  is  - 

a  =  1.53551604307  radii 
e  =  0.00832494892587 
i  =  9 1?14969 13060 
Q  =  224?818733520 
to  =  89?7043902194 

T  =  4356.17523035  days  since  1950.0 
N 

Lj=  225? 44 2972790 

a  =  0.0000429512482776 
xN 

a  =  0.00832483812584 
yN 

|  A_r_|  =  25.42393  kms. 
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h  =  -0.873235460487 
x 

h  =  0.878778312006 

y 

h  =  -0.0248623193590 
z 

t  =  4356.09128442  days  since  1950.0 
N  =  0 

q  =  1.52273295039  radii 
p  =  1.53540962433  radii 
n  =  0.0390833126090  rad/min 
c"  =  0.29  x  10'10  l/mins 


i 


On  the  following  12  pages  the  residuals,  corrections  to  the 

orbital  elements,  and  the  corrected  orbital  elements  for  6  iterations 

are  printed.  The  first  iteration  corrects  n  only  and  the  next  5  iterations 

correct  n,  a  a  „,U  ,fl  ,  and  i. 
xN  yN  o 

The  13th  page  of  the  differential  correction  has  the  following 
information 

Time  of  epoch  =  338.0912842 


Old  rms^  =  8.45  kms 

A  a 

yN 

= 

-0. 28740xl0"9 

_  Q 

New  rms,  =  0.98223  kms 

Au 

= 

-0.31332x10 

1 

Old  rms.  =  0  kms/sec 

o 

An 

s= 

0.7159xl0-1° 

z 

New  rms  =  0  kms/sec 

Ai 

= 

0.1180xl0'10 

An/n  =  0.66858125xl0'9 

Ac" 

= 

blank 

A  a  =  0. 39666x10 
xN 

No.  of  residuals  used  =  13  No.  of  residuals  rejected  =  2 


Old  Elements 
(at  time  of  nodal 
crossing  closest 


to  epoch 


New  Elements 
(at  time  of  nodal 
crossing  closest 


to  epoch) 


New  Elements 


>rediction  rev. 


Rev .  no . 

1288 

Set.  no. 

12 

L(deg . ) 

136. 

TNo(days) 

338. 

a(radii) 

1. 

e 

0. 

i(deg) 

91. 

n(deg) 

224. 

dXdeg) 

89. 

1288 

12 


19780 

136.48960 

06358 

338.06354 

5355160 

1.5395001 

00832 

0.01085 

149 

91.151 

817 

224.813 

735 

93.642 

CQ(days)  0  0 

H  (st.  miles)  2071.6  2071.8 

q 

Pg(min)  160.819  161.445 

=  0.044  deg/day,  CD  =  -1.098  deg/day 


1288 

13 

136.48960 

338.06354 

1.5395001 

0.01085 

91.151 

224.813 

93.642 

0 

2071.8 

161.445 


I 


I 


Two  sets  of  seven  standard  element  cards  with  the  new  elements 
are  punched;  the  first  set  is  for  the  time  of  nodal  crossing  closest  to 
epoch,  and  the  second  set  is  for  the  time  of  nodal  crossing  into  the  pre¬ 
diction  revolution. 

The  first  set  is  as  follows: 

All  cards  have  Satellite  no.  =  163,  element  set  no.  =  12,  and  card  type  —  E 

Card  1  Card  2 

Sat  name  =61  SIGMA  1  Year  of  =  1961 

No 

N0  =  1288  TNq  =  338.0635468638  days  since  1961.0 

e  =  0.0108581954  L  =  136°  .489606287 

o 

i  =  91  .151658969 

Card  4 

c  =  0.000000000029  days 
ftQ  =  0.044228236  deg/day 
<i>o  =  -1.098043234  deg/day 


Card  6  Card  7 

aQ  =  1.5395001036  radii  Nf  =  1288 

Standard  brightness  =  Expiration  date  of  Bulletin  =  1960.0  (meaningless) 

0.000000000 

0^  =  0,0000000000  IST0P  =  0 

The  second  set  is  the  same  as  the  first  set  except  that  the  expiration 

date  of  the  Bulletin  is  December  4,  1961,  l^1  3lm  30s  .45  and  the  element  set 
no.  =  13. 


Card  3 

PNq  =  0. 1121147138  days 

ft  =  224°. 813728127 
o 

O)  =  93°. 642762117 
o 

q  =  1  522783910  radii 
o 

Card  5 

All  blank. 
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CASE  2 


INPUT: 


Control  Card  1 
Sat.  No.  =  034 
Element  set  no.  =  099 
Sat.  name  =  60  EPSI  1 
Epoch  rev.  no.  =  8670 
DC  option  =  0 

Control  Card  3 
Prediction  rev.  no.  =  8670 
i  option  =  0 

Observation  Cards 


Control  Card  2 

c  =  -0.11251x10 

o 

M  =  99 

=  -001  mins 

e.  =  100  kms. 

4 

m  =  2 


Time  of  Obs 


oca. 

No. 

Year 

Mo. 

Day 

Hr. 

Min. 

Sec. 

h(deg) 

A(deg) 

P(kms) 

0316 

1 

11 

25 

4 

10 

26.0 

5.8 

328.75 

1834.67 

0316 

1 

11 

25 

5 

45 

6.0 

5.8 

318.75 

1760.54 

0316 

1 

11 

26 

3 

05 

24.0 

1.6 

321.25 

2149.71 

OUTPUT: 

The  three  position  fixes: 


Time 


YMMDDHH 

MMSS  SSS 

x(radii) 

y( radii) 

z( radii) 

r (radii) 

61112504 

1025.980 

0.0908385425456 

-0.445447206264 

0.963478327438 

1.06534752121 

61112505 

4505.987 

0.267520165478 

-0.415019701875 

0.939419196686 

1.06128074461 

61112603 

0523.986 

-0.0972024607181 

-0.444682810921 

0.958790981211 

1.06135350652 
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Two  possible  solutions  are  printed,  the  first  is  the  true 


solution. 

First  Solution  - 

a  =  1.05957282532  radii 
e  =  0.0227595478820 
i  =  64°. 8369177419 
ft  =  186°. 156960786 
oo'=  197°. 405908140 

Tn  =  4347.22091013  days  since  1950.0 
L  =  280°. 9 74228400 

a  M  =  -0.0217173765006 
xN 

a  =  -0.00680827280506 
yN 

I  Ar  |  =  2506.851  kms. 

Second  Solution  - 

a  =  1.09746945963  radii 
e  =  0.0384727603814 
i  =  115°. 163082256 
ft  =  6°. 15696106455 
oo  =  45°. 7917417283 
Tn  =  4347.26605452  days 

L  =  78°. 6716120957 

a  M  =  0.0268258411204 
xN 

a  „  =  0.0275776638241 
yN 

|  Ar  1  =  5181.064  kms 
100  kms. 


h  =  -0.0998981992772 
x 

h  -  0.926056950702 

y 

h  =  0.437564682332 

z 

t^  =  4347.17391151  days  since  1950.0 

N  =  1 

q  =  1.03545742691  radii 

p  =  1.05902396957  radii 

n  =  0.0681830911664  rad/min. 
c"  =  0. 190794808906x10” 7  1/mins. 


h  -  0.101620034315 
x 

h  =  -0.942018373054 

y 

h  =  -0.445106502215 
z 

t^  =  4347.17391151  days  since  1950.0 
N  =  0 

q  =  1.05524678016  radii 

p  =  1.09584503690  radii 
n  =  0.0646821189089  rad/min. 
c"  =  0. 171704492357x10" 7  1/mins 


since  1950.0 


No  set  of  7  standard  element  cards  should  be  punched  because 
for  both  possible  solutions. 


> 
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SECTION  6 


EXPERIMENTATION 


Tables  I  through  VII  contain  the  results  of  several  cases.  The 
results  also  include  a  direct  comparison  with  the  SPADATS  elements.  Tables 
II,  VI,  and  VII  use  simulated  observations  and  therefore  give  an  excellent 
indication  of  the  true  accuracy  of  the  program.  The  experimentation  contains 
cases  with  a  range  of  eccentricities  from  0.005  to  0.2  and  with  up  to  13 
revolutions.  In  all  cases,  the  number  of  revolutions,  N,  and  the  time,  At, 
between  the  two  observations  is  indicated.  =  0.001  mins,  for  all  cases. 

In  general,  for  each  possible  number  of  revolutions  between  the 
first  two  position  fixes,  there  is  both  a  direct  and  retrograde  two  body 
orbit  that  passes  between  the  two  positions  in  the  proper  duration  of 
time  The  third  vector  check  is  useful  for  isolating  the  true  solution 
orbit.  Comparing  each  possible  residual ,  |Aj:  |  ,  with  a  fixed  is 

not  the  most  efficient  way  to  isolate  the  true  solution  from  the  spurious 
ones.  It  is  nearly  impossible  to  assign  a  value  to  that  will  be 

larger  than  the|A  r|  for  the  true  solution  and  yet  smaller  than  the 
|Ar  s  for  all  the  other  possible  solutions.  For  instance,  if  the  third 
observation  is  from  the  same  part  of  the  orbit  as  the  first  two  (and 
this  is  quite  likely  if  the  third  observation  is  from  the  same  station) 
the  | A  r [ 1 s  for  many  of  the  possible  solutions  can  be  as  small,  or  smaller, 
than  the  |Ar  |  for  the  true  solution.  Similarly,  if  the  third  observation 
is  far  removed  from  the  first  two  in  time,  the  |  Ax  I  can  be  Ruite  larSe 
for  even  the  true  solution.  In  the  bottom  two  cases  of  Table  II,  some 
of  the  spurious  solutions  yielded  a  smaller  |Ar  |  than  did  the  true 
solution. 


A  better  way  to  isolate  the  true  solution  is  to  pick  the  orbit 
with  the  smallest  |A£.|  .  This  does  not  always  yield  the  true  solution, 
as  witness  the  aforementioned  two  cases.  However,  in  all  the  other  cases 
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performed  in  the  experimentation,  the  smallest  |^ir  |  was  associated  with 
the  correct  solution.  Picking  the  orbit  with  the  smallest  l^^jr  |  is 
reliable  if  the  three  observations  are  chosen  judiciously.  The  first  two 
positions  should  be  picked  from  the  same  revolution,  if  possible,  and  if 
not  3  as  close  together  as  feasible.  The  third  observation  should  be 
from  a  different  station  than  the  first  two  in  order  for  it  to  be  on 
a  different  part  of  the  orbit. 

An  examination  of  the  cases  performed  in  the  experimentation 
indicates  that  seme  cases  show  a  rather  large  deviation  from  the  SPADATS 
data  for  the  two  parameters  e  and  co  •  This  is  due  to  the  fact  that  there 
is  a  singularity  in  the  method  for  e  =  0.  Therefore,  if  the  orbit  is 
close  to  being  circular,  there  is  a  likiihood  of  errors  being  introduced 
into  the  results.  These  errors  may  be  aggravated  by  two  parameters.  If 
sin  (v2  -  vx)  is  very  close  to  zero  (i.e.,  v2  -  v 1  <4°,  and/or 
p/r;L  is  close  to  unity  the  errors  will  be  magnified.  For  example,  in 
Table  II  for  "60  Nu  2 ,  A t  -  1.470  mins,"  the  value  for  sin  (v2  -  v^)  is 
quite  small  and  the  error  in  is  37°.  On  the  other  hand,  for  the  same 
satellite  and  At  =  3h  52m  .904,  sin  (v2  -  vj^  is  quite  large  and  the 
error  in  A)  is  1°.  Obviously,  there  is  nothing  that  can  be  done  about 
p/r^  to  improve  the  accuracy,  but  sin  (v2  -  v^)  can  be  improved  in  some 
cases.  If  there  is  a  choice  in  the  selection  of  the  observations,  the 
two  that  give  as  different  positions  on  the  orbit  as  possible  should  be 
chosen.  They  should  both  be  on  the  same  revolution,  if  possible,  in  order 
to  be  consistent  with  the  suggestions  concerning  the  magnitude  of|Ai|  . 

If  this  is  done,  and  p/r^  does  not  turn  out  to  be  very  close  to  unity, 
the  program  will  handle  successfully  orbits  with  very  low  eccentricities. 

The  smallest  error,  for  any  of  the  cases  shown  here,  in  the 
semimajor  axis,  a,  or  perigee  distance,  q,  is  .0004  radii.  To  achieve 
this  same  accuracy  in  p,  the  experimentation  indicates  that  £3  could  be 
as  large  as  0.5  mins.  This  would  save  approximately  10  iterations  as 
compared  with  using  €  3  =  0.001  mins. 
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TABLE  I 


Element 

a  (  r  ad  i  1 ) 
e 

i(degrees) 

.17  (degrees) 

AJ  (degrees) 

T  (days  since 
n  1961.0) 


Element 

a(radii) 

e 

i(degrees) 

J1  (degrees) 
(degrees) 

T  (days  since 
n  1961  0) 


Element 

a( radii) 
e 

i(degrees) 

,1 1  (degrees) 

'  (degrees) 

T  (days  since 
n  1961.0) 


Element 

a( radii  ) 
e 

i( degrees) 

!  I  (degrees) 

A)  (degrees) 

T  (days  since 
n  1961.0) 


60  Epsilon  1 , 

N  =  1,  At  =  94.7  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1.0404 

1.0355 

0 . 0049 

0.01546 

0.02276 

0.00730 

65.018 

64.837 

0.181 

186.424 

186.157 

0.267 

211,092 

197.406 

13.686 

329.22115 

329  22091 

0.00024 
(0 . 34  mins) 

60  Epsilon  1, 

N  =  3,  At  =  284.1  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1  0404 

1.0440 

0 . 0036 

0.01546 

0.01252 

0.00294 

65  018 

65.114 

0,096 

186.645 

186.032 

0.613 

211.120 

195.002 

16.118 

329.15741 

329. 15744 

0.00003 
(0 . 04  mins) 

60  Eta  2,  N  = 

6,  A  t  =  709 . 5  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1.0961 

1.0994 

0.0033 

0.03076 

0.02807 

0.00269 

66.770 

66.687 

0.083 

217.514 

216.714 

0.800 

220. 194 

228. 182 

7.988 

330. 15650 

330. 15654 

0.0004 
(0.06  mins) 

61  Sigma  1,  N 

=  10,  At  =  1739.8  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1.5259 

1.5267 

0.0008 

0.00923 

0.00898 

0  00025 

91- 131 

91.155 

0.024 

223.586 

223.646 

0.060 

106.290 

116.180 

9.890 

315.17846 

315.17836 

0  00010 
(0 . 14  mins) 
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TABLE  II  61  SIGMA  1 


N  =  0,  A  t  = 

10.0  mins 

Magnitude 

Element 

SPADATS 

D- Iteration 

of  Difference 

a(radii) 

1. 5398 

1.5355 

0.0043 

e 

0.01069 

0.00832 

0.00237 

i(degrees) 

91. 198 

91.150 

0.048 

Q (degrees) 

224.778 

224.819 

0.041 

c,  (degrees) 

93.384 

89.704 

3.680 

T  (days  since 
n  1961,0) 

338. 17589 

N  =  1,  At  = 

338.17523 

165.0  mins 

0.00046 
(0.67  mins) 

Magnitude 

Element 

SPADATS 

p- Iteration 

of  Difference 

a( radii) 

1. 5398 

1.5404 

0 . 0006 

e 

0.01069 

0.01152 

0.00083 

i(degrees) 

91.  198 

91.149 

0.049 

Si  (degrees) 

224. 778 

224.979  ‘ 

0.  201 

(degrees) 

93.384 

85.428 

7.956 

T  (days  since 
n  1961.0) 

338. 17589 

N  =  2,  A  t  = 

338.17577 

325  0  mins 

0.00012 
(0.17  mins) 

Magnitude 

Element 

SPADATS 

p- Iteration 

of  Difference 

a(radii) 

1  5398 

1 . 5404 

0.0006 

e 

0.01069 

0.01152 

0.00083 

i(degrees) 

91.  198 

91.148 

0.050 

SI  ( degrees) 

224  778 

225. 232 

0.454 

(degrees) 

93.384 

85.091 

8.293 

T  (days  since 
n  1961  0) 

338  17589 

N  =  4,  A  t  = 

338.17577 

645.0  mins 

0.00012 
(0.17  mins ) 

Magnitude 

Element 

SPADATS 

p- Iteration 

of  Difference 

a(  rad ii) 

1.5398 

1. 5404 

0.0006 

e 

0.01069 

0.01180 

0.00111 

i(degrees) 

91.198 

91.156 

0.042 

'  l  (degrees) 

224.778 

223.458 

0.320 

(degrees) 

93.384 

103.821 

10.437 

T  (days  since 
n  1961,0) 

338  17589 

338,17565 

0.00024 
(0.35  mins) 
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Element 

a(radii) 

e 

i(degrees) 

Q  (degrees) 
l O  (degrees) 

T  (days  since 
n  1961.0) 


Element 

a(radii) 

e 

i(degrees) 

<f?(degrees) 

60  (degrees) 

T  (days  since 
n  1961.0) 


TABLE  II  61  SIGMA  II 
(Continued) 

t  =  1090.0  mins 


SPADATS 

1.5398 
0.01069 
91. 198 
224.778 
93.384 
338.17589 


p- Iteration 

1. 5404 
0.01191 
91.246 
224.826 
95.972 
338.17576 


N  =  8,  At  =  1295.0  mins 


SPADATS 

1.5398 

0.01069 

91.198 

224.778 

93.384 

338.17589 


p- Iteration 

1.5405 

0.01167 

91.146 

225.487 

100.110 

338.17567 


Magnitude 
of  Difference 

0.0006 
0.00122 
0.048 
0.048 
2.588 
0.00013 
(0. 19  mins) 


Magnitude 
of  Difference 

0.007 
0.00098 
0.052 
0.709 
6.726 
0.00022 
(0.32  mins) 
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TABLE  III 


60  Nu  2,  N  =  0,  At  = 

1.470  mins 

Magnitude 

Element 

S PAD ATS 

p- Iteration 

of  Differenc 

a(radii) 

1,1674 

1.1799 

0.0125 

e 

0,01894 

0.01760 

0.00134 

i(degrees ) 

28.22 

28.22 

0.00 

l (degrees ) 

284.883 

284.068 

0.185 

3  (degrees) 

281.228 

318.774 

37.546 

;  (days  since 
n  1962,0) 

5.99040 

5.99182 

0.00142 
(2.04  mins) 

60  Nu  2,  N  -  2,  At 


3h  52.m904 


Magnitude 


Element 

S PAD ATS 

p- Iteration 

of  Difference 

a(radii) 

1.1674 

1.1664 

0.0010 

e 

0.01894 

0.02123 

0.00229 

i(degrees) 

28.22 

28.14 

0.08 

J2  (degrees) 

283.891 

283.542 

0.349 

A,'  (degrees) 

281.215 

282.417 

1.202 

T  (days  since 
n  1962.0) 

5.99040 

60  Xi  1,  N  =  0, 

5.99035 

At  =  1. 133  mins 

0.00005 
(0.07  mins) 

Magnitude 

Element 

SPADATS 

p- Iteration 

of  Difference 

a( radii) 

1.2104 

1.2132 

0.0028 

e 

0. 12000 

0. 11816 

0.00184 

i(degrees) 

49.92 

49.89 

0.03 

( } (degrees ) 

301.317 

301.603 

0.286 

A  (degrees) 

231.657 

226.343 

2.314 

T  (days  since 

26.14511 

25. 14532 

0.00021 

n  1962.0) 

60  Xi  1,  N  =  0, 

At  =  1.231  mins 

(0.31  mins) 

Magnitude 

Element 

SPADATS 

D- Iteration 

of  Difference 

a(radii) 

1.2104 

1.  1970 

0.0134 

e 

0.12000 

0. 13250 

0.01250 

i(degrees ) 

49.92 

50.21 

0.29 

Jl (degrees) 

301. 587 

301.038 

0.549 

(degrees) 

231.432 

234.458 

3.026 

T  (days  since 
n  1962.0) 

26.06703 

26.06543 

0.00160 
(2.30  mins) 
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Element 

a( radii) 
e 

i(degrees) 
iQ  (degrees) 
(degrees) 

T  (days  since 
n  1962.0) 

Element 

a(radii) 

e 

i(degrees) 

■  2  (degrees) 

SO  (degrees) 

T  (days  since 
n  1962.0) 

Element 

a(radii) 

e 

i(degrees) 
(degrees) 
v  (degrees) 

T  (days  since 
n  1962.0) 


Element 

a(radii) 

e 

i(degrees) 
fl  (degrees) 

/,)  (degrees) 

T  (days  since 
n  1962.0) 


TABLE  IV 


60  Pi  1,  N  = 

0,  At  = 

1.122  mins 

Magnitude 

SPADATS 

p-  Iteration 

of  Difference 

1.1058 

1.1060 

0.002 

0 . 00946 

0.00960 

0  00014 

48.57 

48.69 

0.  12 

8.344 

7.911 

0.433 

250.369 

236.382 

13.987 

10.24803 

10.24793 

0.00010 
(0.15  mins) 

60  Pi  1,  N  = 

13,  At  = 

=  22h18I!'876 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1.1058 

1.0876 

0.0182 

0.00947 

0.02407 

0.01460 

48.57 

51.04 

2.47 

327.264 

22.324 

55.060 

287. 285 

223.865 

63.420 

19.11103 

19.11051 

0.00052 
(0. 74  mins) 

60  Iota  1,  N 

=  0,  A  t 

=  1.296  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1.2373 

1.2244 

0.0129 

0.05751 

0.05900 

0.00149 

47.28 

46.74 

0.  54 

356.493 

355.233 

1.260 

271.611 

264.449 

7.162 

30.08049 

=  o,  At 

30.07904 

0.00145 
(2.09  mins) 

60  Iota  1,  N 

=  0.813  mins 

Magnitude 

SPADATS 

p- Iteration 

of  Difference 

1. 2373 

1.2630 

0.0257 

0.05751 

0.02969 

0.02782 

47.28 

47.86 

0.  58 

355.959 

354.653 

1.306 

272.123 

290.813 

18.690 

30. 24188 
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30. 24456 

0.00268 
(3.86  mins) 

Element 

a( radii) 
e 

i(degrees) 
f}  (degrees) 

4'  (degrees) 

T  (days  since 
n  1962.0) 


Element 

a(radii) 

e 

i(degrees) 
(degrees) 
id  (degrees) 

T  (days  since 
n  1962.0) 

Element 

a(radii) 

e 

i(degrees) 

2  (degrees) 

61  (degrees) 

T  (days  since 
11  1962.0) 


TABLE  V 


60  Iota  2,  N  =  0, 

S PAD ATS 

1.2499 
0.01176 
47.22 
35. 263 
167.909 
30. 12975 


At  =  0.888  mins 

p- Iteration 

1.2579 
0.02182 
47.34 
35.647 
319. 155 
30. 13093 


60  Iota  5,  N  =  0,  At  =  0.910  mins 

p- Iteration 


S PAD ATS 

1.2523 
0.01212 
47.20 
47.520 
159. 549 
30. 13276 


1.2891 
0.03581 
47.77 
48.682 
88.026 
30. 13657 


61  Alpha  1,  N«0,At«  1.520  mins 


S PAD ATS 

1 . 0800 
0.00518 
97.40 
113.970 
58.058 
16.21038 


p- Iteration 

1.0900 

0.01527 

97.10 

113.683 

59.571 

16.21128 


Magnitude 
of  Difference 

0.0080 
0.01006 
0. 12 
0.384 
151.246 
0.00118 
(1.70  mins) 


Magnitude 
of  Difference 

0.0368 
0.02369 
0.57 
1.162 
71.523 
0.00381 
(5.48  mins) 


Magnitude 
of  Difference 

0.0100 
0.01007 
0.30 
0.287 
1.539 
0.00090 
(1.29  mins) 
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Element 

a(radli) 

e 

i(degrees ) 

J?  (degrees) 
tl)  (degrees) 

T  (days  since 
n  1962  0) 


Element 

a(radii) 

e 

i( degrees) 

SI  (degrees) 

&  (degrees) 

T  (days  since 
n  1962.0) 


Element 

a(radii) 

e 

i(degrees) 

SI  (degrees) 

£0  (degrees) 

.T  (days  since 
n  1962.0) 


Element 

a( r adii) 

e 

i(degrees) 

' (degrees) 

.  (degrees) 

T  (days  since 
n  1962.0) 


Element 

a(radii) 

e 

i(degrees ) 
l (degrees) 
(degrees) 

T  (days  since 
n  1962.0) 


TABLE  VI  SIMULATED  OBSERVATIONS 
N  =  0,  ^t  =  2.0  mins 
True 

Elements  p-  Ite  ration 


1.5 
0.  2 
30.0 
29.352 
31.028 


1.5046 
0.19717 
30. 189 
28.895 
32.548 
66.43099 


N  =  0,  At 
True 

Elements 

1.5 

0.2 

30.0 

29.352 

31.028 


N  =  0,  At 
True 

Elements 

1.5 
0.2 
30.0 
29 .353 
31.028 


4.0  mins 


p- Iteration 

1.5011 
0. 19921 
29.997 
29.193 
31.606 
66.43081 

10.0  mins 


p- Iteration 

1.4996 

0.20010 

30.023 

29.153 

31.334 

66.43069 


N  =  0,  At  = 

True 

Elements 

1.5 
0.2 
30.0 
29.352 
31  028 


N  =  0,  A  t  = 

True 

Elements 

1.5 
0.  2 
30.0 
29.352 
31.028 


.  0  mins 


p- Iteration 

1.5011 

0.19946 

29.981 

29.218 

31.570 

66.43082 

.0  mins 


p- Iteration 

1.5011 
0. 19936 
29.968 
29 . 225 
31.738 
66.43084 


Magnitude 
of  Difference 

0.0046 
0.00283 
0. 189 
0.457 
1.520 


Magnitude 
of  Difference 


0.0011 

0.00079 

0.003 

0.159 

0.578 


Magnitude 
of  Difference 

0.0004 

0.00010 

0.023 

0.200 

0.306 


Magnitude 
of  Difference 

0.0011 
0.00054 
0.019 
0.  134 
0  542 


Magnitude 
of  Difference 


0.0011 
0  00064 
0.032 
0. 127 
0  710 
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TABLE  VII  SIMULATED  OBSERVATIONS 


N  =  1,  At  =  22.0 

mins 

True 

Magnitude 

Element 

Elements 

p- Iteration 

of  Difference 

a(radii) 

1.5 

1.4992 

0 . 0008 

e 

0.2 

0.19854 

0.00146 

i(degrees) 

30.0 

29.997 

0.003 

(degrees) 

29.676 

29.785 

0.109 

i a)  (degrees) 

30.514 

31.169 

0.655 

T  (days  since 

66.21533 

n  1962.0) 

N  =  2,  At  =  32.0 

min 

True 

Magnitude 

Element 

Elements 

p- Iteration 

of  Difference 

a(radii) 

1.5 

1.4991 

0.0009 

e 

0.2 

0. 19916 

0.00084 

i(degrees) 

30.0 

29.919 

0.081 

fl  (degrees) 

29.676 

29.429 

0.247 

60  (degrees) 

30  514 

31.036 

0.522 

T  (days  since 

66.21527 

n  1962.0) 

N  =  4,  A  t  =  15.0 

mins 

True 

Magnitude 

Element 

Elements 

p- Iteration 

of  Difference 

a(radii) 

1.5 

1.4991 

0 . 0009 

e 

0.2 

0.19988 

0.00012 

i(degrees) 

30.0 

29.217 

0.783 

J l (degrees) 

29 . 838 

29.031 

0.807 

t\\  (degrees) 

30.257 

31.573 

1.316 

T  (days  since 

66.10752 

n  1962.0) 

N  =  10,  At  =  1.0 

mins 

True 

Magnitude 

Element 

Elements 

p- Iteration 

of  Difference 

a( radii) 

1.5 

1.4992 

0.0008 

e 

0.2 

0.19990 

0.00010 

i(degrees) 

30.0 

31.097 

1.097 

^(degrees) 

29.838 

31.313 

1.475 

60 (degrees) 

30.258 

30.435 

0.177 

Tn(days  since 

66.10792 

1962.0) 
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SECTION  7 


PROGRAM  STATUS 


7 . 1  LIMITATIONS 

a.  Inclination  must  not  be  equal  to  zero. 

b.  The  orbit  must  have  an  eccentricity  greater 
than  zero  out  equal  or  less  than  V , 95 , 

c.  The  time  of  the  two  observations,  which  determine 
the  orbit,  must  not  be  equal  to  zero. 

d.  The  difference  between  the  true  anomalies  of 
the  two  observations  must  not  be  equal  to 
zero  or  even  multiples  of  n  . 

e.  The  time  between  the  two  observations  must  be 
less  than  two  days. 

7.2  IMPROVEMENTS  UNDER  CONSIDERATION 

a.  As  was  indicated  in  Section  6,  the  correct 
selection  of  the  three  observations  to  use 
for  this  orbit  determination  scheme  is  impor¬ 
tant  for  correct  solutions.  There  is  no 
reason  why  this  selection  cannot  be  done  inside 
the  program  rather  than  leaving  it  up  to  an 
operator  or  analyst. 
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b.  Taking  the  solution  with  the  minimum  |  A  rj  as  the 
true  solution  would  be  a  better  way  of  isolating 
the  true  solution  than  comparing  the  various 
|At|' s  with  a  fixed  £4. 

c.  Changing  the  equation  for  P2max  to  read: 

P2max  <  1.975  x  (the  smaller  of  r^  or 

would  conform  with  the  second  limitation  mentioned 
above  . 

d.  Max  N  can  be  constantly  updated  to  conform  with 
the  value  of  the  previously  selected  p.  This 
will  reduce  the  number  of  spurious  solutions. 
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